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Abstract 

A careful analysis is made in studying X(4260) property by fitting J/ip7r + TT~ and tt + it~ 
invariant mass spectrum. For the latter final state theorem for tttt, KK couple channel sys- 
tem is used. As a consequence the ratio 9xj/tpa/9xj/ipf (980) * s est i mate d- Meanwhile we 
get a sizable coupling between X and coxco- Our analysis suggests that the X(4260) res- 
onance is a conventional cc state renormalized by LUXcO continuum and the wave function 
renormalization constant is estimated to be Zx — 0.8. The renormalization effect is im- 
portant in reducing the leptonic decay width of X(4260) and to escape from the existing 
experimental constraint. 

PACS numbers: 13.20.Gd, 14.40.Gx, 13.75.Lb. 

Key words: X(4260), charmonium, final state interactions, partial wave analysis 

1 Introduction 

A resonant state X(4260) (previously called Y(4260)) has been found by BABAR collabo- 
ration in initial-state radiation (ISR) process in 2005 [1 ], in J/ip7r + 7i~ final state. The mass and 
width are given as M = 4259 ±8(stat.)tl(sys.)MeV and T = 88±23(stat.)^(sys.)MeV. The 
branching ratio is determined to be 

r e+e - x Br(X -> tt+tT J/tjj) = 5.5 ± 1.0±g;| eV . (1) 

Also the X(4260) state has been confirmed by CLEO [ ] and BELLE [3] experiments, and CLEO 
also has seen it in 7r ir J/ip and K + K~ J/ip processes. BELLE collaboration also gave the 
7r + 7r~ invariant mass spectrum. More recently, BABAR has restudied the J/?pir + ir~ channel and 
obtained upgraded results. [4] 
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Figure 1: X(4260) and nearby resonances from naive quark model calculation [5]. 
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Figure 2: X(4260) and its nearby thresholds. 

In theory aspects, the existence of X(4260) is very interesting, because it is generally 
thought there is not enough unassigned vector states in charmonium spectrum, also the masses 
(including the recently reported Y(4360), X(4630)/Y(4660) states) are inconsistent with naive 
quark model predictions. [ ] The only such 1 states expected up to 4.4 GeV are IS, 2S, ID, 
3S, 2D and 4S, and they seem to be well established. [ ] The situation is depicted in Figure 1. 
It tends to believe that the discrepancy between naive quark model prediction and the observed 
spectrum is ascribed at least partly to the existence of many open charm thresholds, and the lat- 
ter will distort the spectrum. The situation is depicted in Figure 2. However the experimental 
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searches seem to indicate a rather negative result. For example, a candidate final state for X(4260) 
decay is DD*, However, BABAR collaboration could not find X(4260) in DD*, DD*, D*D* [7]. 
BELLE collaboration measured the process e + e~ — > D°D*~n + , and there was no evidence for 
X(4260). [8]. Rather extensive searches in open charm final states in e + e~ annihilation at en- 
ergies between 3.7GeV and 5.0GeV are made, but only upper limits on X(4260) decays are 
given. [9] Furthermore BES collaboration did not find the signal of X(4260) by scanning R value 
as well [10]. 

Because of above reasons many theoretical papers have been produced devoted to the inves- 
tigation on X(4260), and it is generally believed that the existence of X(4260) signals a degree 
of freedom beyond conventional cc state. In the literature, many proposals have been made, 
e.g., charmonium [1 1], XcoP° molecule [12], uxd molecule [13], ccg hybrid state [14, 15], A C A C 
bayronium [16], D\D or D Dq molecule [17, 18], etc.. 

In this article, we also devote to the study of X(4260) issue, inspired by the recent BABAR 
experiment. [4] Through a careful analysis to experimental data available, we find that the 
X(4260) has a sizable coupling to cuXco channel, but not to other (nearby) channels. The paper is 
organized as the following: Section 1 is the introduction, in section 2 we review on theoretical 
tools we use in this paper, special emphasis is made on the final state interactions between pions. 
In section 3 we give a detailed description to our numerical fit program. In section 4 we extract 
the X coupling to J/ipa and J ftp / (980) and determine the pole location of X(4260) on Rie- 
mann sheets. In section 5 we estimate the wave function renormalization constant of X(4260) 
induced by the interaction between X,lu and Xco- Finally in section 6 we draw conclusions based 
on the studies of this paper. 

2 Theoretical Discussions on X — > J/i/jtttt Process 

2.1 Effective Lagrangian Describing X — > J/ ipirir Interactions 

Assuming that the A (4260) state being a J PC = 1 chiral singlet, we write down the 
following effective lagrangian describing X — > J/ipmv process, which is accurate in the leading 
order in the expansion in terms of n momentum in the center of mass frame of nn system, 

Cx^pp = hxX^ < u a u a > +h 2 X^ < X+ > +h 3 X IMX ^ < u p u a > , (2) 

where we have used the anti- symmetric representation (X^, ip^ u ) to describe the 1 state 
X(4260) and J/ip, respectively, and F^ u denotes the photon field strength. Notice that in the 
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present notation, one has, 



e+e ~ 3 M x ' K } 



when neglecting electron mass. Also notice that, up to 0(p 2 ) level, in Eq. (2) there exist only 
three independent interaction terms with coefficients hi, h 2 and h 3 . Further, = i(u + d^u — 



ud fl u + ) and 



is a parametrization of the pseudo-goldstone octet: 

/>° + > - + K+ \ 



(5) 



The chiral symmetry breaking term with coefficient h 2 in Eq. (2) reads, 

X+ = u + xu + + ux + u, x = 2B diag(m u , m d , m s ) . (6) 

Parameters F n and B can be fixed phenomenologically: F w ~ 92.4MeV and < | | >= 
— F 2 B [1 + 0(m q )}. The Eq. (2) can also be rewritten in an explicit form, 

£i = ^X^ u F^(d p 7i + d p n~ + ^d p n°d p 7r° + d p K + d p K~ + d p K°d p K° + \d p i)dPf]) (7) 



£ 2 = -^X^iml'K+K- + -m^rV + m 2 K K+K- + m 2 K K°K° + {Un\ - -m%) m ) (8) 



+ -dpK°d a K° + -dpK°d a K° + -dprfd^rf) . (9) 
2 2 2 

It is noticed that, according to Eq. (2), only s-wave and rf-wave of the tttt system are involved 
(more discussions will be given later). Also Eq. (2) does not take into account of possible strong 
final state interactions between two pions in the decay product. More discussion on this topic 
will be left in section 2.3. 



2.2 Kinematics and tree level amplitudes 

Denote the momentums of e~, e + , X(4260), J/ip, tt + and 7r _ as qi, q 2 , q, qo, q + and q~ 
respectively, as seen from figure 3. The polarization of J/ip is represented as e A , and define 
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Figure 3: A depiction of kinematics. 

k± = q + ± g~, then one has the following relations, 

s = k 2 ,, 



k_ = —sp = Am — s, (p 

V s 

ql = M 2 j^ k + -k_=0, 
1 

2*,y - ivi j/4> 



k + -q = -(q 2 -M 2 -s). (10) 



The amplitude of X — > J /ipix + n process at tree level is 

*A tree = M Itln («2^ {QU shxu(q 2 , s>) {[Ah.hs - 2ml) + ^m\\{q° ■ qe x - q ■ e^q 



+\h 3 [ -q° a q ■ e 4 ,{k x + k% - k x _k a _) + e x q° a q^k + - k a _k P _) 

-q°\ a e^{k a + ki - k a _kt) + q° ■ qe^ a {k x + k% - k x k a _) ] }, (11) 

where Dx(q 2 ) is the denominator of the X(4260) propagator. In the simplest form, it can be 
parameterized as D x {q 2 ) = M\ — q 2 — iM x T X - 

Follow the helicity amplitude decomposition method, [19] choose the basis of tensors 

t (0) = 1, (12) 
= F, (13) 
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fPO = - = <T - , (14) 

3 k\ 

it is then easy to recognize that terms proportional to hi and h 2 in Eq. (11) only contribute to the 
s-wave whereas the term proportional to h 3 contributes both to s-wave and <i-wave. The overall 
s-wave amplitude reads, 

iAf ee = 111" r M qi, shxu(q 2 , s')e^ { [^-(s - 2m 2 ) + ^ 2 ml](q° • q - q 6 q ox ) 

Mj/^q 2 D x {q 2 ) 2 

+\h,[-\p 2 {s)qlq e s + (l - lp 2 ( S )) kUV ■ k + + \p 2 {s)q« ■ q g™ s 
~ (l " lp 2 (s)) 9»K -qk + .q°- \p\sWs + (l - ^ 2 (*)) k\q° e q ■ k + 
+ \p 2 {s)q° ■ q g xe s -(l- ±p 2 (s)^ k\k%q ■ q° } }, (15) 
whereas the <i-wave part is 

- a tree ^^^^?0 — / \ / f\ 

lAd = Mj^qW^f^ S)7AM(92 ' S ^ 

\h 3 (+t Xa q° a q e - tfg xe q% + tfg pe q x q a - t Xa g ae q° ■ 

~ F^Mj/^q 2 D x {q 2 ) ' ^ 

The square of d-wave amplitude is simple (if neglecting final state interaction effects) after inte- 
grating over angle #* defined in Figure 3, 

f 1 dcosBl-Y,\B d \ 2 (17) 
J— i „/ 



^pl[M%^s - V) + Mj^(6g 4 - q 2 s + s 2 ) - M% (4g 6 + qU - 34gV + 9«* 
+Ml + (q 2 -s) 2 (q 4 + 3q 2 s + 6s 2 )] 



Since the <i-wave contribution is proportional to the 4-th power of the kinematic factor, it is highly 
suppressed comparing with the s-wave contribution. It will be shown later through numerical 
studies that the <i-wave contribution is roughly less than 1% of the total decay rate. 

2.3 Final state interactions 

The tree level amplitude as described in section 2.2 is not sufficient to describe the X — > 
J/ipinr decay process, since the nit final state is mainly in the IJ=00 system and undergoes strong 
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final state interactions (FSI). To include FSI we propose the following decay amplitude: [21] 

Aj =AT e a 1 {s)T ll {s) + Al ree a 2 (s)T 21 (s) , 

A 2 = AT e a l {s)T l2 {s)+^ e a 2 {s)T 22 {s) , (18) 

where subscripts 1,2 denote nix and KK final states, respectively. Especially T n , T 12 and T 22 
represent ixix — > nn, nn — > KK, KK — > KK scattering amplitudes, respectively. Functions a; 
are mild polynomials which place the role to offset the 'left hand' cuts on the complex s plane 
in amplitude T that would not appear in function A. Expressions in Eq. (18) are remarkable in 
the sense that the unitarity relations that the decay amplitudes A- L - as an analytic function of s - 
have to obey: 

ImAi = AXpxTu + A* 2 p 2 T 21 , 

lmA 2 = A* lPl T 12 + A* 2 p 2 T 22 (19) 

are automatically satisfied. Noticeably in a\(s) an additional pole term is added: 

(i) 

ai (s) = — ^ — + + 4 1} s + ■ ■ ■ , (20) 

s- S A 

where s A represents the Adler zero of Tn. The role of the pole term is to cancel the Adler zero 
hidden in T n but not welcome in A. [21] An advantage of the pole term in Eq. (20) is that, by 

(1) c (1) 

appropriately choosing coefficient c as lim s ^ SA Tn(s) = 1, it guarantees 

A x = AX ee + 0(s 2 ) (21) 
in the small s region, which allows the standard chiral power counting to work here. 

2.4 The 7T7T scattering T matrix 

We choose three different txtx scattering T matrices to perform three different fits to the 
X(4260) spectrum and the invariant mass spectrum of the txtx final state: couple channel Pade 
approximation (named as fit I, parameters listed in Table 1), [20] K-matrix unitarization (named 
as fit II, K 3 parametrization of Ref. [21], using parameters refitted in Ref. [22]), and the PKU 
representation (fit III). [23] For more details about these T matrices, we refer to the original 
literature. Choosing three different T matrices is not redundant. The PKU representation gives 
a reliable description to the low energy physics related to the a meson, but lacks of information 
on full couple channel dynamics and hence physics related to / (980). The other two methods, 
on the other side, are not accurate when discussing the a meson but have full respects to couple 
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Li 


L 2 


L 3 


U 


U 


2L 6 +L 8 


10.32 


13.30 


-38.03 


3.770 


17.15 


10.39+4.96 



Table 1 : Low energy constants from couple channel Pade amplitudes of fit I. Here these param- 
eters are refitted and are slightly different from Ref. [20]. The unit is 10~ 4 . 



Pade K-matrix PKU 

sheet-II 0.459 - 0.229i 0.549 - 0.230i 0.463 - 0.247i 

sheet-II 0.989 - 0.013* 0.999 - 0.02H 0.980 - 0.019* 

sheet-Ill - 0.705 - 0.327* 

sheet-Ill - 0.977 - 0.060*' 

Table 2: Poles' location on the -y/i-plane, in units of GeV. 

channel unitarity and hence may provide a better description to the physics related to /o(980). 
Pole positions of the three different T matrices are listed in Table 2. Regarding to the Adler 
zero discussed in section 2.3, In this work, we use sa and cjj as calculated from each T matrix 
(listed in Table 3), instead of using the ^PT 0(p 2 ) value to fix the parameters %\— 167nF 2 ) and 
sa(= "4/2) in the I,J=0,0 channel. 2 





Pade 


K-matrix 


PKU 


Alder zero s 

c (1) 


0.00954 - 0.00015i 
0.32956 -O.OOllOi 


0.01880 
0.14132 


0.00383 
0.36088 



Table 3: The Alder zero location on s-plane and Cq (in units of GeV 2 ). 

3 Numerical Studies and Discussions 

3.1 Introduction to the fit program and nearby cut singularities 

The X(4260) may have other decay channels besides the J/ipini(KK) one. It may also 
couple to a channel that it does not decay for kinematical reason. The nearby thresholds are 
X c0 w (4197MeV), D* + D*~ (4224MeV), 7J-/J 1 + (2420) (4291MeV), X a u (4293MeV), etc.. It 

2 Allowing these two parameters to vary does not make any significant changes on other physical output in the 
numerical fit discussed in section 3. 



8 



is possible that X(4260) couples to these channels. To be specific we attempt to write down the 
denominator of the propagator of X(4260) as, 

D x (q 2 ) = M 2 x -q 2 - iy/tf{ gi k x + g 2 k 2 + T(q 2 ) + T ), (22) 

where k x (k 2 ) is the cm. momentum of channel 1 (2) that lies below (above) the peak position 
around 4260MeV. For example, channel 1 may denote D* + D*~ or coXco, channel 2 may denote 
D~ Df (2420) or Xci w. The latter two are indistinguishable in our study since they are both s- 
wave decays and the thresholds are very close to each other. T(q 2 ) is calculable (see for example 
Eq. (17)) which includes both J/ipirir and J/ipKK decays, T is a constant and simulates the 
decay width into other (un-observed) channels. 

In our program we fit the X(4260) line shape in the region 4.15GeV< < 4.47GeV, 
using data from Ref. [3] (17 data points) and Ref. [4] (17 data points) (see Figure 4a). If one 
tries to fit a larger q 2 , one encounters the problem that s becomes too large and hence invalidates 
the parametrization of low energy T matrices. 3 For txtx invariant mass spectrum we use the 
data given in figure 4b of Ref. [3], which corresponds to m^+^-j/^ e [4.2,4.4]GeV (17 data 
points), and the data £ [4.15, 4.45]GeV (41 data points) from Ref. [4]. In Ref. [3], a MC study of 
efficiency correction at q 2 = 4.26GeV is also given. Through a numerical test we however find 
that the efficiency curve is well reproduced by the txtx two body phase space up to a normalization 
constant, hence in our fit we simply use the two body phase space instead of the efficiency 
corrected one. We also assume the BABAR data of Ref. [ ] have a similar behavior. There is 
no KK invariant mass spectrum available. But one has the integrated J /ijjKK events and this 
contributes another data point, a(K + K~ J/ip) = ± lpb. [2] Summing up there are totally 
93 data points to be used. 

Parameters needed in our fit are the following: Firstly, Eq. (2) describing tree level X(4260) 
J /if) txtx(KK) interactions and the 7 - X transition provides 4 parameters. Secondly, it is found 
that taking 011,2(5) to be linear polynomials (except the Adler zero term) is already good enough 
for data fitting, hence the two oti(s) contribute another 4 parameters. 4 Thirdly, through rather ex- 
tensive numerical tests in different environment, we find that the coupling constant g 2 in Eq. (22) 
is always vanishingly small, therefore we often freeze parameter g 2 in the fit. Furthermore, in 
the absence of experimental evidence on X(4260) decay into other channels we often assume T 
being vanishing. 5 So we often switch off these two parameters unless stated otherwise. Hence 

3 For \fcp- = 4.47GeV, it corresponds to a value sfs ~ 1.17GeV which is within the range that T can provide a 
reasonable description. 

4 In PKU parametrization there is no a-2 so there are only two parameters here, only tree level amplitude of 
X -> J/ ipKK are used, 

5 Adding To will actually not lead to any important change in physical conclusions in this paper. 
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Eq. (22) contains another 2 parameters. Finally, there are two normalization factors N 1 and N 2 
for Belle and Babar ix + ix~ invariant mass spectrum. So summing up there are totally 12 (10 for 
PKU representation) parameters. 

To begin with, we assume channel 1 represents D*D* and channel 2 represents DDi(2420) 
(and/or cu Xci)- In this case it is noticed that the X — > D*D* decay is p wave, hence in Eq. (22) 
ki should be replaced by k\. It is found that all three fits give a fair description to the data. The 
parameters are determined and listed in Table 4 (w/o T ) and Table 5 (with r ). In the fit we 
find that, as already mentioned before, the coupling constant g 2 is vanishingly small. It suggests 
that the coupling of X(4260) to DDi(2A20) and coxa is negligible. On the other side, the 



Fit I (Pade) 



Fit II (K-matrix) 



Fit III (PKU) 



Xd.o.f 

2o(MeV) 
9i 

92 

My(GeV) 



125.5 
93-13 

2.582 ±0.105 
0.76xl0~ 8 ± 0.0005 
0.69xl0~ 7 ± 0.00001 

4.261 ±0.005 



127.8 
93-13 

2.616± 0.090 
0.54xl0- 3 ±0.003 
0.43xl0- 3 ±0.0003 

4.261 ±0.005 



150.2 
93-11 

2.944 ±0.218 
0.17xl0~ 6 ± 0.0005 
0.57xl0~ 7 ± 0.0001 

4.269 ±0.007 



Table 4: Fit results assuming an X(4260) coupling to D*D* S and DDi, T not included in the fit. 



Fit I (Pade) 



Fit II (K-matrix) 



Fit III (PKU) 



Xd.o.f 

So(MeV) 
9i 

92 

My(GeV) 
r (GeV) 



125.0 
93-14 

3.358 ±0.199 
0.14xl0~ 6 ± 0.001 
0.18xl0~ 8 ±0.0002 
4.249±0.004 
0.098 ± 0.010 



127.3 
93-14 

2.836± 0.085 
0.74xl0~ 6 ±0.0006 
0.80xl0~ 7 ±0.0001 
4.257±0.006 
0.020± 0.018 



149.9 
93-12 

3.745 ±2.088 
0.14xl0" 6 ± 0.0006 
0. 13 xl0~ 6 ± 0.00002 
4.260±0.013 
0.052± 0.077 



Table 5: Fit results assuming an X(4260) coupling to D*D* and DD\\ T Q included in the fit. 



coupling gi is also negligible, but the reason is different from neglecting g 2 . Unlike the s-wave, 
the XD*D* coupling is p wave, and this contribution in Eq. (22) is severely suppressed. 
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3.2 Coupling to ujXcQ channel 

Remember that the to Xco channel is in s-wave, we therefore redo the above analysis by 
letting channel 1 being tu Xco- 6 The result is plotted in Figure 4 and Table 6. Comparing Table 6 
with Table 4 (5) we find that including coupling the coXco channel significantly improves the fit 
quality. Meanwhile, the pole position of X(4260) is determined on sheet III of yfs plane, i.e., the 





fit I (Pade) 


fit II (K-matrix) 


fit III (PKU) 


Xd.o.f 

So(MeV) 


108.3 
93-12 

33.416 ± 19.828 


107.8 
93-12 

34.404 ± 19.391 


135.2 
93-10 

34.592±14.944 


9i 


0.476±0.056 


0.446 ±0.048 


0.507± 0.060 


My(GeV) 


4.274 ±0.007 


4.271 ±0.005 


4.279±0.007 




54.031 ±4.958 


55.326±4.647 


52.494 ± 4.858 


N 2 (pb- 1 ) 


19.103 ±1.775 


19.761±1.771 


16.909±1.672 


h 


0.0542 ±0.033 


0.107±0.051 


0.01 8± 0.009 




-0.554±0.496 


-1.077±0.520 


-0.085±0.037 




-0.004± 0.048 


0.0007 ±0.025 


-0.014±0.008 


cf) 


0.633± 0.430 


0.321 ±0.212 


5.876 ± 1.908 


c? 


-0.271 ±0.282 


-0.155 ±0.110 


-3.341 ±1.155 


c (2) 


-0.038±0.011 


-0.003 ±0.014 





c (2) 

c 2 


0.029±0.013 


-0.0006±0.013 






Table 6: Parameters given by fit to X(4260) data, assuming an XluXco coupling g-y. 

one seen by experiments. The results are listed in Table 7. It is remarkable to see that the (real 
part of the) pole is very close to (or even below) the coXco threshold, and be much lower than the 
peak position. There is another pole, lying roughly in similar position, on sheet IV. On the other 
side, from Figure 4a) one finds the BABAR data peak at 4.227GeV. The discrepancy between 
the peak position and the pole position is due to the complicated analytic structure around the 
uXco threshold. To further test the importance of the ujXcqX coupling gi, we also provide the fit 
result by switching off g\ in Table 8 and in Figure 5. The pole position when eliminating g x is 
also looked for and the result is given in Table 9. 

At this stage let us firstly summarize what we have obtained: We find the coupling to chan- 
nel tuXco very important in influencing the fit quality and the output physical picture related 

6 For the parameter Tq and 52 in Eq. (22), it is noticed that including them does not improve the fit quality and 
does not alter any of the major conclusions given below. Also the fit in this subsection gives vanishingly small 
values of Tq and g-2- For this reason, we no longer discuss these parameter hereafter. 
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E(Gev) ITI (Gev) 




(c) 

Figure 4: A fit to the X(4260) invariant mass spectrum by assuming that it couples to cu Xco- 



to X(4260). When turning on the coupling g v , one gets a large g as well, and hence a large 
T x ^e+e-- That means the Br(X — > J/tpim) is very small. As a result, poles on sheet III and 
sheet IV are very close to each other numerically and both are very close to the coxco threshold at 
4197MeV. The picture may suggest that X(4260) be a u>XdO molecular resonance. If one however 
switch off the ojXco channel, one gets a quite different scenario: Coupling g becomes very small 
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Fit I 


Fit II 


Fit III 


sheet-Ill 


4236.6-63.4i 


4237.5-58.9i 


4236.1-68.5i 


sheet-IV 


4237.0-63.0i 


4237.8-58.6i 


4236.6-68.0i 



Table 7: The pole location of X(4260). Assuming an XuXco coupling g 1 . 



as seen from Table 8, which corresponds to T e + e - ~ lOeV. Meanwhile X decays dominantly 
into J/ipirn and the pole locates at around 4250MeV. Confidentially, the x 2 fit favors the former 
scenario. Comparing Table 6 with Table 8, it is found that the total x 2 increases rather rapidly, 
which corresponds to a significance level for parameter g 1 being 3.6a (I), 4.3<r (II) and 4.2o" (III), 
respectively. This strongly suggests that X(4260) has a sizable coupling to channel coXco- 1 

However, before jumping to a quick conclusion we notice that there is a potentially severe 
problem regarding to the results given in Table 6. The central value of go given there, which 
characterizes the transition between X(4260) and the photon field, or T X -> e +e- according to 
Eq. (3), is too large to avoid the upper bound established by an analysis on the R value measured 
by BES experiment: Ref. [24] gives T x ^ e +e- < 420eV at 90% confidence level. Nevertheless 
the value of g given in Table 6 contains a large error bar as well. It is actually understood well 
why such a phenomenon occurs. When g being large, parameters hi have to become small to 
accommodate for the experimental value given in Eq. (1), as can be understood from Eq. (15). 
This results in a non-important T(q 2 ) contribution to the denominator of X(4260) propagator 
(see Eq. (22)). As a consequence, x 2 becomes inert with respect to the change of go when it is 
large enough since the effect can be counterbalanced by a change of hi. The situation is clearly 
depicted in Figure 6 which exhibits an approximate scaling law for parameter g . This means 
that we can not reliably determine parameter g at all when g is large enough, though the value 

7 If one writes down the effective lagrangian of X(4260) coupling to loXcO in the anti-symmetric form 

C = giY^co^Xco , 

then if taking gi ~ 0.5 according to Table 6 we get an estimation on g\ ~ 7.2GeV. 





fit I (Pade) 


fit II (K-matrix) 


fit III (PKU) 


Xd.o.f 


121.6 
93-11 


126.2 
93-11 


152.9 
93-9 


So(MeV) 


2.529 ± 0.087 


2.533 ±0.095 


2.790 ± 0.168 


My(GeV) 


4.260 ±0.005 


4.259 ±0.005 


4.264±0.005 



Table 8: Parameters given by fit to X(4260) data with g x switched off. 
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(c) 

Figure 5: A fit to the X(4260) invariant mass spectrum with parameter g\ switched off. 

of g x is very stable when varying g in Figure 6. It is important to notice from Figure 6 that there 
exists a large enough space for g to maintain a (almost) minimal x 2 while being easily escapable 
of the BES T x ^e+e- bound given in Ref. [ ]. As an example, we list the fit result in Table 10 
by fixing g at 9.62MeV which corresponds to T x ^ e +e- — 211eV. The pole location in such a 
situation is listed in Table 1 1 . 
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Fit I 


Fit II 


Fit III 


sheet-II 


4250.2-51.3i 


4253.0-33.0i 


4256.6-49.1i 



Table 9: The pole location of X(4260) with g x switched off. 




Figure 6: The approximate scaling law of parameter g . 



Before ending this section it maybe worthwhile to briefly comment on the rf-wave contri- 
bution. It is controlled by parameter h 3 in Eq. (11). According to Table 6, h 3 is very small and 
probably be consistent with 0. The estimate is however very unstable, maybe partly due to the 
fact that we do not yet take into account d wave final state interactions. [25] It is also worth 
mentioning that our analysis reveals that there ought to be non-negligible KK — >■ txtx (and vice 
versa) final state interaction. This can be understood from Table 4 by comparing the \ 2 difference 
between the PKU representation and the other two approaches. 

At last we would like to point out that, we also tested the possibility whether the above fit 
results be distorted by possible interfering background. To do this we add a background term by 
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changing the X propagator: 



1 



M 2 x -q 2 -t^( gi k 1 + T(q 2 )) 

1 a e i(t> 



M\-q 2 - iy/tf{g x k x + T(g 2 )) M 2 - q 2 - iM x Y } 



+ c 2 . (23) 



Taking fit I for example, in doing this we observe that adding 6 additional parameters (c 2 is 
complex) only reduces total x 2 by 3. Moreover, parameter g x only decreases very slightly. The 
only important change comes from g - its central value is now roughly 4.1MeV. Hence we 
conclude that the qualitative picture obtained in previous discussions of this section is rather 
robust. 





fit I (Pade) 


fit II (K-matrix) 


fit III (PKU) 


Xd.o.f 


108.9 


108.3 


136.3 


93-11 


93-11 


93-9 


3o(MeV) 


9.615 


9.615 


9.615 


9i 


0.435±0.048 


0.412 ±0.043 


0.450± 0.048 


My(GeV) 


4.273 ±0.006 


4.270±0.005 


4.277 ±0.007 


Ndpb- 1 ) 


54.239± 4.954 


55.511±5.140 


52.799± 4.953 


N 2 (pb- 1 ) 


19.147 ± 1.807 


19.795±1.935 


16.951±2.053 


h 


0.187 ±0.063 


0.383± 0.255 


0.063± 0.065 




-2.030±0.798 


-3.835±2.126 


-0.296±0.149 


h 3 


0.010± 0.151 


0.002±0.0773 


-0.049±0.007 




0.583± 0.477 


0.313 ±0.411 


5.928 ±7.26 




-0.250± 0.271 


-0.150± 0.195 


-3.360 ± 4.354 


(2) 

c\ 


-0.031 ±0.055 


-0.003 ±0.022 





c (2) 


0.022±0.056 


-0.0006±0.021 






Table 10: Parameters given by fit to X(4260) data, assuming an XuXco coupling g\ and keeping 
g fixed. 
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Fit I 


Fit II 


Fit III 


sheet-Ill 


4239.1-61.2i 


4239.9-57.2i 


4240.9-66. li 


sheet-IV 


4243.3-56.1 i 


4243.2-52.9i 


4246.5-59.0i 



Table 11: The pole location of X(4260). Assuming an XuXco coupling gi, g fixed at 9.62MeV. 
The results of this table is similar to Table 7. 



4 Couplings gx^a and g X m /o ( 9 80) 



The analytic continuation of amplitude A - as a function of complex variable s - to different 
sheets can be obtained 8 : 

-2ipiT/ 2 



A 11 = AB n = A ( 1 + 2 ^ T i / i i+2*pit£ 



( l+2ip 2 r 2 J 2 -2i P iT[ 2 
dct S ActS 
dct S det S 



A IV =AB IV = A( _ 2i ^ r i ° 

\l+2ip 2 T' 2 / 2 l+2ip 2 T 2 / 2/ 

where ^4 = (Al, *A 2 ), and T matrix expressions on different sheets are similar, 

T u — TBn , T IH = TBin , T IV = TB IV ■ (25) 

With these expressions the residues or coupling constants of second and third sheet poles can be 
extracted. In the vicinity of pole position on sheet II, one has 

r" (7r7r _> 7r7r)(s _> ^) = Tll( * fj) , (26) 

^(^jj) (a - 2//) 

A n (X _> 2// ) = ^ l( * Jj) (27) 

<S n (2//)(s - ^//) 



and through the definition of coupling constants: 



g 2 

T(im — > mr)(s — > zr) = — ; 

Z R ~ S 

A(X -> ^tttt) (5 -»• z fl ) = g JSUllUL (28) 

- s 



For more details on how to obtain these formulas one is referred to Ref. [26]. 



17 



one obtains the expressions of gf n7T and gxvj for second sheet pole f: 



9 



Tu(z n ) A^zn) 



9fnn ~ ~s[J^d ' 9xvf9f ™ ~ ~s[J^7) ' (29) 

Denote the third sheet pole as z in , which is the zero of det S, then 

det S(s ->■ 2m) = (det S)'(> m )0 - z m ) . (30) 

Coupling constants are hence obtainable, according to Eq. (25) and Eq. (28): 

2 S 22 (z HI ) 2 Sn(zin) -Tnjzni) m , 

^ " 2i Pl dfitS'(z m ) ' ~ 2ip 2 detS'(z III ) ' " detS'(z m ) ' W 

From Eq. (24) one further gets 

Ai{z ni )S 22 {z HI ) A 2 (z III )2ip 2 (z III )Ti 1 (z III ) 

9x *t X9 *™- (detS)W) (detS)'(*m) 

A(z//j)2z>i(zj/j)T 1 J 2 (z m ) A2(zin)S 11 (z In ) 

gXVf X fl^tf - — - — r — , (32) 

(det S)'{z IH ) [det Sy(zin) 

or more easily, by making use of Eq. (18), 

9xvf = AX ee {z n Mx{zui)9f^ + A 2 ree (z in )a 2 (zni)gfKK ■ (33) 

These expressions enable us to get the same value of gx^f, f° r me third sheet pole /, as expected. 

From the method as mentioned above, we can predict the X(4260)-J/ 1 I / -/o(600) coupling and 
the X(4260)-J/\E r -/o(980) coupling. However, owing to the approximate scaling law discussed in 
section 3.2, the magnitude of gxvf is not physical and only the ratio gx^a/gxwf (980) makes 
sense in the scaling region. Table 12 provides an example of such calculation. It is seen however 
that this ratio is not quite stable and is rather sensitive to different fit method. Nevertheless one 
expects that one can evaluate the order of magnitude of the ratio through Table 12. For example 
both Fit II and Fit III give a value of l^^l a few times larger than |<7xw/o(980)l' nence me 
X(4260) — > J/ipirir decay mainly goes through the a or / (600) intermediate state. This agrees 
qualitatively with the estimation [4], 

Br(X (4260) -» J/^/ (980), / (980) -> tt+tt-) 

Er(X(4260) ^ J/1> " ' 



9 One can instead get expressions for g 2 f KK and gx^ / x g/KK and deduces the same expression for gxy f, as it 
should. 
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Pade 


K-matrix 


PKU 




sheet-II 


sheet-II 


sheet-Ill 


sheet-II 


q 2 

JCTTVTT 


-0.122 -0.1 27i 


-0.0727 -0.1 69i 




-0. 147-0. 133i 


n 2 


-0.024 -0.016i 


-0.063 -0.071 i 






9f (980)irn 


-0.043 -0.002i 


-0.071 -0.011 i 


-0.100+0.022 i 


-0.041 -0.005 i 


9 2 fo(980)KK 


0.151-0.026 i 


-0.101 +0.085 i 


-0.022 -0.095 i 




9XJ/i>a 


-10.102-3.311i 


-6.311-4.045i 




-9.394-3.639i 


9XJ/i.>fo(980) 


-2.845-1. 287i 


-4.398-0.607i 


1.383+4.425i 


-6.676-1.233i 


9XJ/4,f () (9S0) 


2.510+2.299i 


1.533+0.708i 




1.458+0.276i 



Table 12: The coupling of X(4260) to / (980) and a (or / (600)). The last three columns provide 
the results by taking g = 9.62. 



5 An Estimation of Wave Function Renormalization Constant 
of X(4260) 

Coupling constant g is related to the two point proper vertex describing 7-X transition. 
After considering the interaction between X, Xco and u, g gets renormalized, 

<7* = Z)?gS , (35) 

where g$ is the value of go at tree level, which may be recognized as the value obtained from 
simple potential model calculation without considering the continuum mixing, and g$ is the one 
measured experimentally in the leptonic decay process. On the other side the wave function 
renormalization constant can be expressed as 



'X 



(36) 



l-Re£'(/z 2 ) ' 

where S is the self-energy, /j, is the physical (pole) mass. According to the approximation made 
in previous sections, 

S(g 2 ) ~ -i^gth . (37) 
More accurately it is recasted in an expression using dispersion relation, 

SO,.) = I f (38) 
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from which one further deduces 



Im E(g' 2 ) 

J(M x +m») a (<l' 2 ~ V 2 ) 



(39) 



The above integral is divergent in relativistic theory and hence not calculable. However, in the 
non-relativistic approximation, it becomes calculable for s wave interaction: 



Im£(g 2 



^Qih -)> -| v^M^^g 2 - (M x + mJ2) 



(40) 



and bringing it back to Eq. (39) one gets 



Zx 1 



1 -ReS^/i 2 ) = 1 



9i 



Re 



V 7 ^ 



(41) 



where M f 2 ft = (M x + m u ) 2 . If X is a bound state then 

1 1 



'X 



l-Re£V) i + ^^nf 



(42) 



637MeV. If X is a 



where we have let [i = M th — e and e is the binding energy, m R - 
resonance then when fi ~ M th — iT/2 one gets instead of Eq. (42) the following formula, 



1 



j x 



l-Re£V) l + 



(43) 



If for example taking gi ~ 0.435 and the third sheet pole of Fit I from Tables 10,1 1, it is estimated 
-^x(4260) — 0.74 using Eq. (43), and Z x — 0.82 using Eq. (41). The potential model calculation 
of go does not take into account the hadron loop effect and hence only corresponds to a 'tree level' 
calculation. The uxco loop correction leads to that the 'tree level' value of T e + e - be reduced by 
a factor Z x ■ 

Finally we comment a few words on X(3872) before ending the present discussion. Con- 
sidering the first row of Table 3 in Ref. [27] for example. From which one gets g\ = 0.0937, 
and here m R ~ 966.5MeV, T ~ 0.4MeV. Bring them back to Eq. (43) one gets Z X (3872) — 0.38, 
which is rather close to the one obtained in Ref. [28]. In obtaining this number other contribu- 
tions like J/ipp and J/ipui are not yet considered, hence one concludes that the mixing effects 
from continuum in X(3872) is stronger than in X(4260). Axiomatic quantum field theory indi- 
cates that < Z < 1, and Z = for a pure composite particle. In this sense, one may conclude 
that X(4260) is more elementary than X(3872). 
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6 Conclusions 



To conclude, in the analysis made in this paper we find a sizable Xu Xco coupling. Other 
couplings to nearby thresholds are found vanishingly small with no significance importance. For 
example, there is no evidence for a gxD*D* ■ This does not support the proposal of X(4260) being 
a D*D* molecule. We also find no evidence to include an XDD\ coupling. This is a rather 
negative observation towards the hybrid assignment of X(4260), since the latter indicates a large 
hybrid coupling to a spin and spin 1 meson. [29] The large coupling to Xco w is also of spin 0, 1 
type, but a ccg hybrid should couple much stronger to a charmed meson pair than to a cc state 
plus a light meson, since the latter needs a color rearrangement. Also, the observation made in 
this paper seems to exclude the DD\ (coxa) assignment of X(4260). Finally the nix spectrum is 
apparently and overwhelmingly in isoscalar s wave, this certainly excludes the proposal that X 
be a XcoP molecule. 

Of course once the g\ coupling is switched off the X(4260) behaves very differently: It has 
only a tiny electron decay width but with a dominant J/ipim decay mode. Nevertheless, this 
picture is disfavored by the x 2 analysis. After all, the numerical study of this paper suggests 
that the X(4260) state contains a sizable coXco component, and one may attempt to describe 
X(4260) as an coXco molecule, or more appropriately, a molecular resonance. Nevertheless, this 
seemingly natural conclusion is unlikely for two reasons. Firstly the coupling g\ is still far 
from being large enough to regard X(4260) as a pure molecule. The second evidence against the 
molecule interpretation comes from the value of g given in Table 6 and Figure 6, where it is seen 
that the fit prefers a larger g , even though we can not actually determine g reliably. A larger 
go easily leads to a T e + e - of a few hundreds eV, which is hard to understand in a pure molecule 
picture. It is actually highly plausible that X(4260) contains a large portion of cc component. The 
preference for a large go, in our understanding, seems to support the cc origin of X(4260). In 
Ref. [30] a screening potential is used to calculate the charmonium spectrum and it is found that 
the mass of 4 3 Si state is around 4273MeV with a r 4 3 Si~>e+e- — 970eV. The latter is dangerously 
large considering the BES bound [24]. Taking into account S-D mixing effect may roughly 
reduce the above leptonic width by half [ ] . Here it is pointed out that the renormalization from 
tuXco continuum can further reduce T e + e - of the bare 4 3 S' 1 state roughly by a factor 0.8. This is 
helpful in allowing X(4260) as a candidate cc state to escape from the constraint set up by BES 
experiment [24]. Furthermore, a mixture with coXco continuum is certainly helpful in reducing 
the decay rates into open charm channels as well. 
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